Is there spin-charge separation 
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The spin and density correlation functions of the two-dimensional Hubbard 
model at low electronic density (n) are calculated in the ground state by using 
the power method, and at finite temperatures by using the quantum Monte 
Carlo technique. Both approaches produce similar results, which are in close 
agreement with numerical and high temperature expansion results for the two- 
dimensional t — J model. Using perturbative approximations, we show that the 
examination of the density correlation function alone is not enough to support 
recent claims in the literature that suggested spin and charge separation in the 
low electronic density regime of the t — J model. 
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The normal state of the high temperature superconductors does not behave 
as an ordinary Fermi liquid (FL). 1 For the last several years, Anderson 2 has 
strongly supported the idea that instead it may be described as a Tomonaga- 
Luttinger liquid (TLL). 3 The TLL state has been shown by exact theories 4 and 
numerical studies 5 ' 6 to be the ground state of the Hubbard model, and also 
the t — J model, in one dimension. A special property of the TLL state is the 
separation of spin and charge degrees of freedom. 

For these ideas to be applicable to the cuprates, the crucial question is 
whether phenomenologically realistic two dimensional models of correlated elec- 
trons present features in the ground state similar to those of their one dimen- 
sional counterparts. Recently, several numerical studies have addressed this 
important issue, namely the possibility of spin-charge separation in the two- 
dimensional t — J model. Using a ground-state projection technique ("power" 
method) to study the low electronic density region, two of us 7 have found indi- 
cations that spin, charge and pairing correlations behave in a similar manner as 
in ID. Qualitatively the TLL state seems to provide a consistent phenomeno- 
logical interpretation of the numerical data. In parallel, based on the results of 
high temperature expansions, Putikka et al. 8 have argued that spin-charge sep- 
aration occurs in the 2D t — J model at low and high electronic densities (par- 
ticularly at high electronic density) . The gauge theory 9 approach also predicts 
a non-Fermi liquid (NFL) behavior, although it may not necessarily correspond 
to a TLL state. On the other hand, for the 2D Hubbard model, analytical 
studies 10 based upon diagrammatic methods suggest the presence of a Fermi 
liquid at low electronic density. 

Since it is well known that the t — J model is equivalent to the Hubbard 
model in the strong coupling limit, the results described above are apparently 
inconsistent. However, the nonperturbative constraint of no double occupancy 
in the t — J model may produce subtle differences with the Hubbard model. 
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Precisely, one of the purposes of this paper is to study numerically the possible 
variation of physical properties between these two models as the constraint of 
no double occupancy is relaxed. Our analysis shows that the density and spin 
correlation functions in the ground state of the 2D Hubbard and t — J mod- 
els are qualitatively similar at least at low electronic density. To examine the 
question of spin-charge separation we compare our results for the density corre- 
lations obtained by the power method at zero temperature on 8x8 and 16x16 
clusters, and by the Quantum Monte Carlo (QMC) method 11 at finite tem- 
perature, with that of the high temperature expansions. 8 Excellent agreement 
between the results of these three methods is obtained. However, from our anal- 
ysis we cannot identify the spinless fermions (SF) Fermi wavevector 2kfi F as 
the characteristic wavevector of the t — J and Hubbard models at low electronic 
density as suggested by Putikka et al. 8 Both the density and spin correlations 
can be understood qualitatively in terms of perturbative approaches such as a 
random-phase approximation (RPA) 12 with a renormalized Hubbard coupling 
U. In addition, the correlations in real space at low electronic density are shown 
to decay so rapidly with distance that the subtle issue of spin-charge separation 
in these models is difficult to address unless these small correlations are accu- 
rately evaluated on a finite cluster, or the functional form of the correlations at 
large distance is obtained with some reliable technique. 

The Hamiltonian for the 2D t — J model considered here has the form 

H t j = -t ]T (c+c j(T + h.c.) + J ^ (Si • Sj - inmj), (1) 

<i,j><7 <i,j> 

with the constraint of no double occupancy. The Hamiltonian for the 2D Hub- 
bard model is well-known and will not be reproduced here. For very large U/t 
the Hubbard model is equivalent to the t — J model with J = 4t 2 /U, up to 
two-particle hopping terms. The good agreement between results for the t — J 
and Hubbard models reported below justifies the omission of these three-site 
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terms. Two numerical methods are used to calculate the equal-time density 
and spin correlation functions, N(q) and S(q), defined by the relations 

N(q) = ^e^ r <5n 5n r >, (2) 

r 

S(q)=^e^ r <S z S z r >, (3) 

r 

where S 2 . = \ Y, a p c+ Q a^ /3 c r/3 , and 5n r = Y.a c ta c va ~ (n). Here (n) is the 
average density of electrons. The brackets in Eqs.(2) and (3) refer to thermal 
averaging in the grand canonical ensemble when the QMC method is used. 
At zero temperature, the ground-state wave function obtained by the power 
method in the canonical ensemble is used to calculate the average. We have 
observed that the well known fermion determinantal sign problem does not pose 
a difficulty in the low-electronic-density region considered in this paper in any 
of the techniques. However, this problem becomes more severe with increasing 
density, and thus we restrict our analysis to the low density region. 

The power method has been proven to be very effective in calculating 
ground-state correlation functions of the t — J model in ID 6 and 2D. 7 The 
ground state wave function is projected out by applying a large power of the 
Hamiltonian, (— H) p , to a trial wave function. The power p required to reach 
convergence depends on the choice of the trial functions. For the case of the 
Hubbard model we use the well-known Gutzwiller wave function, 14 i.e. \GW) = 
g D \FG), where \FG) is the ideal Fermi gas wave function on a lattice and D 
is the total number of doubly occupied sites. Here, g is the only variational 
parameter. At g = 0, the factor g D becomes the well known projection operator 
Pd that projects out states with doubly occupied sites. For the t — J model we 
use the wave function first proposed by Hellberg and Mele 15 in ID and later 
generalized by Valenti and Gros 16 to 2D. This function, which we shall call 
HMVG, is basically of the same form as \GW(g = 0)) i.e. a Slater determinant 
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for up-spin electrons and one for down-spin electrons. In addition to these 
two determinants, it contains a long range correlation part between all the 
particles, n i<J |r« — Yj\ u (while for nearest-neighbor particles we chose v = 0). 
It was shown in Ref.7 that this wave function is very close to the ground state 
for J/t < 2 in the low-electronic-density region. 

S(q) and N(q) are plotted in Figs. la and lb, respectively, against momenta 
along the T-X-M-r directions for an 8x8 cluster. The open circles represent 
the QMC results obtained at temperature T = t/10, U = 4t and (n) = 0.159. 
The open squares are for U = 8t and (n) = 0.155, at T = t/8. The U = 8t 
results deviate further from the ideal Fermi gas shown by the dashed line than 
the U = 4t data. The open triangles represent ground-state results for U = 8t 
obtained for the same lattice with 10 particles ((n) = 10/64 0.156), by 
applying p = 12 powers of the Hubbard Hamiltonian to the |GW) wave function 
with g = 0.5. In most regions of q-space the open triangles and squares take 
the same values. It is gratifying to find out this excellent agreement between 
two very different numerical techniques, namely the power method and QMC. 
To gauge the effect of the constraint of no double occupancy, in the same figures 
we also present results for the t — J model 7 at J = O.lt. These results, which 
are represented by the solid triangles in Fig.la-b, are obtained from the trial 
wave function HMVG-z/ = 0.1 with power p = 16. The t — J model results agree 
very well with the open triangles which corresponds to the Hubbard model with 
U = 8t (J = t/2 in the t — J model language). Hence, there is little difference 
in the correlation functions between strong and intermediate couplings. 

Let us now analyze the implications of Fig.la-b. All the curves in Fig. la 
have peaks at the 2kp wavevectors, except the ideal Fermi gas which only 
presents a discontinuous derivative. The peak size increases with the value of 
U/t. The presence of these peaks 7 implies a stronger spin- density- wave cor- 
relation at finite coupling than the ideal gas. It is interesting to notice that 
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a similar peak is observed in ID. Though the magnitude of the peak in ID 
is much greater than in 2D, this large difference may be partly due to the 
dimensionality effect. Unlike ID in which there is only one 2kp wavevector, 
in 2D there is a characteristic vector in each direction in the two dimensional 
momentum space, each one carrying a peak in S(q). In contrast to the spin 
correlations, N(q) shown in Fig. lb is reduced at q = 2kp as compared to the 
values of the ideal Fermi gas shown by the dashed line. The reduction is larger 
as the coupling U/t is increased. The plateau observed for k > 2k p for the ideal 
gas seems to have shifted to a larger value of k. Putikka et al. 8 have argued 
this new wavevector to be 2k F F . The presence of the Fermi wavevector kf F 
for a spinless fermionic system would imply the separation of charge and spin 
degrees of freedom. 

To examine the important issue of whether 2k F F appears in the numerical 
results, and in order to reduce possible finite size effects we have calculated 
N(q) on a 16x16 lattice using both the QMC and power-method techniques. 
In Fig.2, N(q) is plotted as a function of q along the diagonal direction q x — q y . 
The open squares represent the QMC results obtained at temperature T = t/4, 
U = 8t, and (n) = 0.2. The open triangles represent ground-state results ob- 
tained for the same cluster with 50 particles (i.e. (n) ps 0.195) by applying 
twelve powers of the t — J Hamiltonian to the HMVG variational wave function 
with v = 0.04. The results of the high temperature expansions 8 corresponding 
to J = t/2 and temperature T = J/2 are indicated by the solid line. There is 
very little difference between the results of these three numerical methods. This 
implies that the numerical accuracy of the data is not questionable, but only 
the interpretation needs to be analyzed carefully. For comparison, the result of 
the ideal Fermi gas is plotted as the thin continuous line in Fig.2. It is obvious 
that the results of the interactive system deviates appreciably from the ideal 
gas results. However, these results do not unambiguously support the identifi- 
cation of a singularity at 2kf, F . A more conservative interpretation is that they 
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indicate a broad maximum of N(q) at q = (tt, it) which may just reflect the 
short-range effective repulsion between particles. 

To explore further this assumption we compare N(q) against the T = 
RPA results (dashed line in Fig. 2). The best fit is obtained by choosing the 
renormalized interaction U to be 12t. Thus, the apparent shift of the character- 
istic wavevector can be mimic very well by a simple perturbative (Fermi liquid 
based) approach. We have also calculated S(q) using RPA 12 . The peaks at 2k F 
are also reproduced this time using a smaller effective coupling U = 3t. Such 
a qualitative description of spin and density correlation functions in terms of 
a simple RPA tends to support the point of view that the ground state of the 
Hubbard model is just a strongly correlated Fermi liquid. But this may also be 
misleading. It is also possible to reproduce the correlation functions of a one 
dimensional Hubbard model by using the RPA approximation. Just like in two 
dimensions, a small effective interaction U is enough to produce a large peak 
in S(q) at the proper wavevector 2kp. To fit the density correlation function, 
a larger U is needed i.e. the systematic behavior is very similar in ID and 2D 
(as emphasized in Ref.7). 

The excellent agreement between techniques that work at zero and finite 
temperature shown in Fig. 2 suggests that the shift in N(q) cannot be due to 
subtle long distance correlation functions but to short distance effects. To study 
this hypothesis we analyzed in real space the density-density correlation, C(r), 
for the case of the one band Hubbard model. Fig. 3a shows that this correlation 
decays rapidly with distance and it becomes negligible at four lattice spacings 
away from the origin (numerically the signal at this distance is approximately 
5 x 10~ 4 C(r = 0)). These 2D correlations are considerably smaller than those 
obtained in the case of the one dimensional Hubbard model, which we know 
shows spin-charge separation. This analysis shows that it would be difficult 
to obtain reliable numerical information about the behavior of the correlation 

7 



functions at distances larger than a few lattice spacings. Thus, a proper study 
of spin-charge separation seems beyond present day accuracy of computational 
and series expansion analysis at low electronic density. 

Can the results of our analysis be extended to higher densities? In Fig. 3b, 
N(q) is shown at quarter-filling using the Hubbard model with U/t = 8 and 
the QMC technique. The results deviate considerably from the non-interacting 
Fermi gas, but they can be accurately reproduced by a simple perturbative 
calculation (first order) with an effective coupling U = 4t (see also Ref.17). 
Then, we believe that our conclusions for the Hubbard model can be extended 
to the domain 0.0 < (n) < 0.5. At higher densities the perturbative approach 
breaks down at intermediate and large couplings, due to antiferromagnetism. 
On the other hand, the results for the t — J model at this density are very 
similar to those of non-interacting spinless fermions, as remarked in Ref.8 using 
the high temperature expansions. The possible origin of this discrepancy is 
currently under study. 

In summary, we have presented spin and density correlation functions for 
the one band Hubbard model at low electronic density. The ground-state re- 
sults obtained by the power method agree well with the finite-temperature 
results obtained by QMC. Measurements at intermediate U/t couplings for the 
Hubbard model are consistent with strong coupling data for the t — J model. 
Compared against the ideal gas results, we confirm that N(q) is appreciably 
reduced at q = 2kp as claimed by Putikka et al. 8 This difference increases with 
the strength of the Hubbard interaction. The enhancement of spin-density- wave 
correlation as shown by the appearance of peaks at the 2kp wavevectors for S(q) 
also increases with U/t. This result, first observed in Ref.7, is confirmed by the 
present study on larger clusters and thus finite size effects seem small. On the 
other hand, the RPA approximation can provide a rough qualitative under- 
standing of all these results. In addition, examining N(q) on a 16x16 lattice 
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we did not find evidence for the presence of the characteristic wavevector of a 
spinless Fermion model. Actually, the density correlations in real space decay 
so rapidly that making any statement about their asymptotic behavior based on 
numerical techniques at finite temperature is risky. Thus, based on the current 
available information it is not possible to conclude that spin-charge separation 
takes place in the low electronic density of the 2D Hubbard and t — J models. 
However, we cannot rule out this possibility either. The complete separation of 
spin and charge as in the infinite U limit of the ID Hubbard model may not be a 
proper guidance for 2D studies. A possible scenario is that although charge and 
spin are separated, they interact strongly as in the finite U/t Hubbard model 
in ID. This question is currently being studied. 
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Figure Captions: 

Fig. 1 (a) Spin correlation function S(q) and (b) density correlation function N(q) 
in momentum space along the T-X-M-r directions on a 8 x 8 square lattice. 
Open circles represent QMC results at U = 4t, (n) = 0.159 and T = t/10. 
The open squares correspond to U = 8t, (n) = 0.155 and T = t/8. The 
open triangles represent ground state results for the Hubbard model at 
U = 8t with 10 electrons. The solid triangles denote ground-state results 
for the t — J model at J = O.lt. The dashed line is the result for an ideal 
Fermi gas. 

Fig. 2 The density correlation N(q) along the diagonal direction q x = q y on a 
16 x 16 cluster. The open squares represent the QMC results obtained at 
T = t/4, U = 8t and (n) = 0.2 (results on 12 x 12 lattices are also shown). 
The open triangles are power method ground-state results for the t — J 
model at J = t/2. The solid line is the result obtained by the high tem- 
perature expansion of Ref.8. The dashed line denotes the RPA prediction. 
The result for an ideal Fermi gas is plotted as a thin continuous line. 

Fig. 3 (a) Density-density correlation C(r) = (n n r ) — (n) 2 as a function of the 
distance r. The solid line is the power method result obtained at J/t = 0.5 
(i.e. U/t = 8 in the Hubbard model) on a 16x16 cluster and density 
(n) w .0.20. The dashed line denotes the result for a tight-binding non- 
interacting system U/t = of the same size. The correlations are consid- 
ered along the diagonal of the lattice (the lattice spacing is equal to one); 
(b) The density correlation N(q) along the direction T-X-M-r using the 
QMC technique (full squares) on a 16 x 16 cluster, at T = t/4, U/t = 8, 
and density (n) = 0.5. The dashed line indicates the result for an ideal 
Fermi gas, while the dotted line corresponds to a perturbative calculation 
using an effective coupling U = 4t. The continuous line corresponds to 
non-interacting spinless fermions at the same density and temperature. 
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